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Summary of Research 


In this grant, we have proposed a three-year research effort focused on de- 
veloping High Performance Computation and Communication (HPCC) metnoa- 
ologies for structural analysis on parallel processors and clusters of workstations, 
with emphasis on reducing the structural design cycle time. Besides consolidating 
and further improving the FETI solver technology to address plate and shell struc- 
tures we have proposed to tackle the following design related issues: (a) parade 
coupling and assembly of independently designed and analyzed three-dimensional 
substructures with non-matching interfaces, (b) fast and smart parallel re-analysis 
of a given structure after it has undergone design modifications, (c) parallel eval- 
uation of sensitivity operators (derivatives) for design optimization, and (d) last 
parallel analysis of mildly nonlinear structures. 

While our proposal was accepted, support was provided only for one year. 
During that year, we addressed only the following two topics. 


An optimal FETI solver for realistic aerospace structure s 

Most if not all realistic aerospace structures consist essentially of plate and 
shell components, for which the original (one-level) FETI method is not numeri- 
cally scalable. For this reason, we have first developed a two-level FETI method 
for solving iteratively on massively parallel processors large-scale systems of equa- 
tions arising from the finite element discretization of static and dynamic plate 
bending problems. This method is essentially an extension of the FETI domain 
decomposition algorithm to fourth-order problems. The main idea is to enforce 
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exactly the continuity of the transverse displacement field at the substructure cor- 
ners throughout the preconditioned conjugate projected gradient iterations. This 
results in a two-level FETI substructuring method where the condition number of 
the preconditioned interface problem does not grow with the number of substruc- 
tures, and grows at most polylogarithmically with the number of elements per 
substructure. We have demonstrated these theoretically proven optimal conver- 
gence properties of the new FETI method for several finite element plate bending 
static and transient problems. Then, we have extended this approach to shell 
problems, revisited the preconditioning problem, addressed the efficient imple- 
mentation of the corresponding iterative solvers on massively parallel processors, 
highlighted the computational price of mathematical optimality and discussed its 
consequences, and reported some impressive performance results on two Paragon 
XP/S and IBM SP2 massively parallel processors for several realistic plate and 
shell based aerospace structures. 

Parallel coupling and assembly of three-dimensional heterogeneous 
substructures with non-matching interfaces 

The need for assembling independent finite element substructure solutions 
arises in several aerospace engineering problems including the design and analysis 
of complex structural systems, component mode synthesis, global/local analysis, 
adaptive refinement, and parallel processing. We have addressed the solution of 
such problems by a two-field hybrid method where the substructures are joined 
with low-order polynomial or piece-wise polynomial Lagrange multipliers, and 
present a Rayleigh-Ritz based smoothing procedure for improving the accuracy 
of the computed coupled solution in the presence of various substructure hetero- 
geneities. We have considered both conforming and non-conforming substructure 
meshes, and demonstrated the benefits of the proposed two-step solution method 
with several examples from aerospace structural engineering. 
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